Deformation in ^^Si* produced via ^^O + ^^C reaction 
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The energy spectra of the a particles emitted in the reactions ^^0 (7-10 MeV/nucleon) + ^^C 
have been measured in the center of mass angular range of 25° < 6c.,n. ^ 70°. The experimental 
energy spectra have been compared with those obtained from the statistical model calculation 
with "deformability" parameters predicted by rotating liquid drop model (RLDM) and also fitted 
the same with optimized "deformability" parameters, which are quite different from the respective 
RLDM values. The data have also been found to be explained quite well using "frozen" deformation 
approximation, where the "deformability" parameters have been kept fixed at RLDM values of 
the parent nucleus throughout the decay process. The effective radius in the latter case is smaller 
compared to that obtained using the optimized parameters; however, in both cases, the deformations 
(effective radii) are larger than the corresponding RLDM values. So, considering the uncertainties 
in the estimation of actual compound nucleus deformations, it can, only qualitatively, be said that 
equilibrium orbiting, which is similar to particle evaporation in time scale, could also be one of the 
contributing factors for the observed deformation. 

PACS numbers: 24.60. Dr, 25.70.Jj, 25.70.Gh, 27.30.+t 
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I. INTRODUCTION 

Exploring the role of clustering in fragment emission 
mechanism is a subject of current interest in low energy 
nuclear reaction. Several attempts have recently been 
made to understand the effect of cluster structure on the 
reaction mechanism of light a-cluster nuclei, e.g., ^''Ne + 
[IHi, '^Mg + I, 28Si + 0,0 etc. In each 
case, an enhancement in the yield and/or resonance-like 
excitation function in a few outgoing channels not very 
different from the entrance channel has been observed, 
which was indicative of significant contribution from the 
deep-inelastic orbiting (DIO) mechanism. Through this 
process, a long-lived dinuclear composite is formed 8], 
which is not fully equilibrated in all (in particular, shape) 
degrees of freedom and therefore decays preferentially 
around the entrance channel; at higher excitation energy, 
however, the decay is expected to be of fully equilibrated 
in nature. It is thus interesting to probe into the limit 
of survival of such dinuclear composites at higher excita- 
tions. 

The study of deformation of the excited composite is 
important to differentiate the DIO composite from the 
fully equilibrated compound nucleus (CN) because the 
deformation should be larger in the former case. Detailed 
studies of the deformation of ^"Ca* [9i,|ia] and ^^S* ^ 
were done recently with the help of light charged particle 
(LCP) spectroscopy to extract "anomalous" deformation 
in these cases confirming the presence of DIO mecha- 
nism. For these systems, orbiting ^^Ca* [1, Q, ^^S* 
[ll-Q)had already been conjectured from the fragment 



emission studies. Recently, we have reported that, for 
a-cluster system ^^O + ^^C [l^, there was enhancement 
in the boron product yield at excitation energies ~67-85 
MeV (beam energy in the range of 7-10 MeV/nucleon), 
which indicated the survival of long-lived dinuclear or- 
biting composites at such high excitation energies. An 
orbiting dinuclear complex is assumed to be more de- 
formed than a fully shape equilibrated CN; this moti- 
vated us to look for such enhanced deformation in a ^^O 
-f i^C dinuclear system at same excitation energies us- 
ing the LCP spectroscopy technique. In this technique, 
LCPs, evaporated from the excited composites, are used 
to study the properties of hot rotating nuclear systems as 
a function of excitation energy, angular momentum, and 
deformatio n [ill [l"3l - [25j u sing statistical model codes e.g., 
CASCADE '26^yGANES [IJ, CACARIZO ^]. 

The article has been arranged as follows. The ex- 
perimental setup has been described in Sec. HT] The ex- 
perimental results and analysis have been presented in 
Sec. mil The details of the statistical model calculation 
has been given in Sec. HV] The extracted values of the 
deformation of ^^Si* and their possible implications have 
been discussed in Sec. [V] Finally, the summary and con- 
clusion have been given in Sec. I VII 



* |skundu@vecc.gov.in 

t Present address: Dept. of Physics and Centre for Astroparticle 
Physics and Space Science, Bose Institute, Block EN, Sector V, 
Bidhan Nagar, Kolkata - 700 091, India 



II. EXPERIMENTAL DETAILS 

The experiment was performed at the Variable Energy 
Cyclotron Centre, Kolkata, using ^^O ion beams of en- 
ergies 117, 125, 145 and 160 MeV. Self-supporting 
of thickness ~ 514 /ig/cm^ was used as a target. The a 
particles were detected using a Si(SB) telescope lO/im 
AE, ^ 5 mm E). The typical solid angle covered by the 
telescope was '^0.3 msr. The calibration of the telescope 
was done using elastically scattered ^^O ion from a Au 
target, as well as a ^^^Th-a source. Inclusive energy dis- 
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FIG. 1. Energy spectra (cm.) of a particles obtained at dif- 
ferent angles for different beam energies. The symbols repre- 
sent experimental data. The dash-dot-dashed and solid lines 
represent CASCADE calculations with RLDM and optimized 
values of spin dependent 'deformability' parameters, respec- 
tively (see Table HI]). All experimental and calculated spectra, 
starting from the lowest angle, were multiplied by 10"^ 10"\ 
10°, 10\ 10^, 10^, respectively. 
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FIG. 2. Comparison of the energy spectra (cm.) of a- 
particles obtained at cm. angles 26° (open triangles) and 66° 
(inverted triangles) for beam energy 145 MeV. 



B. Angular distribution 

The cm. angular distributions, {da/dO) cm., obtained 
by integration of the cm. energy distributions for the 
beam energies of 117, 125, 145 and 160 MeV, are sho-wn 
in Fig. [3] as a function of 9c.m.- In all cases, the values 
of {da/d9)c.m. are found to be constant over the -whole 
range of observed cm. angles. So, {da/dfl)c.m. is oc 
l/sin^c.m., which is characteristic of the emission from 
an equilibrated composite nucleus. 



C. Average velocity 



tributions for the LCPs -were measured in the laboratory 
angular range of 11° - 29°, -which covered the angular 
range ~ 25° - 70° in the center-of-mass frame. 



III. EXPERIMENTAL RESULTS 

A. Energy spectra 

The center-of-mass (cm.) energy spectra of the a par- 
ticles at different beam energies are sho-wn in Fig. [TJ Be- 
cause of the inclusive nature of the spectra, there may 
be some admixture of contributions of other direct reac- 
tion mechanism (such as pre-equilibrium emission, etc.) 
to the equilibrium emission spectra at for-ward angles in 
particular. In Fig. [21 the energy spectra obtained at lab 
angles 11° (^c.m.= 26°, open triangles) and 29° (6'c.m.= 
66°, solid inverted triangles) at £^io&=145 MeV are com- 
pared; it is seen that the slopes of the both spectra are 
completely matched -with each other. It is thus clear that 
the effects of other direct reaction mechanisms are not 
significant at the measured angles. 



The average velocities, Vav, of the a-particles emitted 
at different angles have been extracted from the average 
energy, Eav , -which has been calculated using the expres- 
sion 

■where the index i covers the "whole energy spectrum in 
the laboratory. The parallel (w||) and perpendicular {v\^ 
components (-with respect to the beam direction) of Vav 
at each angle are plotted in Fig. |4l It is seen that they 
fall on a circle -with the center at CN velocity, vcn, and 
radius of average velocity in cm., faJ"', "which implies 
that the average velocities (as -well as energies) of the a 
particles are independent of the cm. emission angles. It 
again indicates that the a particles are emitted from a 
fully energy equilibrated source moving "with the velocity 

VCN- 

IV. STATISTICAL MODEL CALCULATIONS 

It is confirmed from the experimental results that the 
Qf-particles are emitted from a fully energy equilibrated 
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FIG. 3. Angular distribution of a particles as a function of 
cm. angles, 6c. m.- Symbols represents the experimental data 
and solid lines show fit to the data obtained using {da /d6)c.m. 
= constant. 



composite which may be a CN. So the statistical model 
code CASCADE 2^ has been used to explain the present 
experimental data because this code is based on the as- 
sumption that the evaporation of particles takes place 
from an excited CN which is in full equilibrium with re- 
spect to all degrees of freedom. 



A. CASCADE calculations procedure 

The spin distribution in CASCADE is usually derived 
from the known fusion cross section using the strong- 
absorption model wherever data are available; otherwise, 
it is calculated theoretically. The decay of the CN by par- 
ticle emission is calculated in the framework of Hauser- 
Feshbach formalism [2^ The probability that the 

parent nucleus (excitation energy Ei , spin Ji , parity tti ) 
emits a particle x with an orbital angular momentum 
kinetic energy ex ( €x^Ei-E2 ), and spin s is given by 

Pxdex = {^xj — r — -, r 

h 27rfi,pi , Ji , TTi ) 



J2+S Jl+S 

X Tiiex)dex, 

S=\J2-s\ l=\Ji-S\ 



(2) 



where E2 , J2 and 1^2 are the excitation energy, spin, and 
parity of the daughter nucleus, respectively, pi and p2 
are the level-densities of the excited parent nucleus and 
daughter nucleus, respectively, and, S = J2 + s is the 
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FIG. 4. Average velocity curve. Solid circles correspond to 
experimental data and solid lines show fit to the data obtained 
using the equation v\ = (va^')'^ - (fn - vcnY where 
is the average velocity of a-particles in the cm. frame. The 
arrow indicates the position of vcn- 



channel spin. The transmission coefficients Ti{ex) for the 
scattering of particle x on the daughter nucleus (inverse 
process) are obtained by using standard optical model 
(OM) potential for elastic scattering. In the OM trans- 
mission coefhcient calculation, the parameters have been 
taken from the references given in the Table HI 



B. Input parameters 

The standard form of CASCADE is quite successful 
in explaining the LCP evaporation in the light ion in- 
duced reaction in general, where the CN is assumed to be 
nearly spherical. However, in case of heavy-ion-induced 
reaction, there is appreciable deviation between the ex- 
perimental and the predicted LCP evaporation spectra. 
This deviation is attributed to the deformation of the 
excited compound system which is angular momentum 
dependent. Therefore, to explain the LCP energy spec- 
tra, the effects of the deformation of the CN should be 
included in the statistical model calculations. The defor- 
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TABLE I. Input parameters used for CASCADE calculations 
for the reactions ^''O + ^^C at beam energies 117, 125, 145, 
and 160 MeV. 



Angular momentum distribution; 

Critical angular momentum Jcr as in Table Ull 
Diffuseness parameter AJ = Ih 

OM potentials of the emitted LCP and neutrons 

(1) Neutrons: Wilmore and Hodgson [311] . 

(2) Protons; Pcrey and Perey [33 |. 

(3) a-particles: Huizenga and Igo [s^ . 

(4) Multiplication factor of the OM radius; REACT = 1 

Level-density parameters at low excitation; E < 7.5 MeV 

(1) Fermi-gas level-density formula with empirical pa- 
rameters from Dilg et al. [3J|. 

Level-density parameters at high excitation; E > 15MeV 

(1) Fermi-gas level-density formula with parameters from 
LDM (Myers and Swiatecki [35| ) 

(2) Level-density parameter a = A/S MeV^^ 

Yrast line 

Igff — Jo(H-5i J*), 5i, S2 are given in Table Ull 

7-ray width (Weisskopf units) 

(1) El = 0.001 

(2) Ml = 0.01 

(3) E2 = 5.0 



mation affects the particle spectra in two ways. First, it 
lowers the effective emission barrier, and second, it in- 
creases the moment of inertia. The first effect modifies 
the transmission coefficients for the evaporated particles 
which may be taken care of by increasing the radius pa- 
rameter of OM potential. However, the change in mo- 
ment of inertia affects the level density and the slope 
of the particle spectrum. This can be taken care of by 
incorporating the spin-dependent deformability parame- 
ters m, Is^. For level density calculations, excited 
energies have been divided into three regions. 

Region I (low excitation energy, i? < 3 to 4 MeV): 
Here, the experimentally known discrete levels are used 
for all the nuclei produced in the cascade. In some cases, 
known high-spin states at higher excitation energy are 
included as yrast levels in region H. 

Region II (medium excitation energy, 4 MeV < E < 
7.5 MeV) : Analytical level density formula is used in this 
region. The parameters a and A are deduced em piri cally 
for each nucleus from the work of Vonach et al. [36| and 
Dilg et al. (s^]. The excitation energy is corrected for 



the parity effects. 

Region HI (high excitation energy, E > Eldm)- Shell 
effects and parity corrections are neglected in this region. 
The same formula is then used but with LDM parameters 
taken from Ref. [35]. 

Between regions II and HI, the level density parame- 
ters are interpolated linearly. The parameters are given 
in Table H] The level density used in regions II and HI 
for a given angular momentum J and excitation energy 
E is given by well-known Fermi-gas expression [TT| with 
equidistant single-particle levels and a constant level den- 
sity parameter a: 



3/2 



1 



2I,ffJ (E + T-A-Ej)^ 
xexp \2{a{E - A - Ej)y^^] , (3) 

where T is the thermodynamic temperature and A is 
the pairing correction. The rotational energy, Ej, is ex- 
pressed as. 



Ej= (7^]JiJ + l)- 



12 J, 



The effective moment of inertia, leff, is written as. 



(4) 



(5) 



where Iq {— ^A^^^Vq) is the rigid-body moment of iner- 
tia. Si and $2 are deformability parameters, tq is radius 
parameter, and a is the level density parameter. So, from 
the above equations, it is clear that by changing tq, a, 
Si, and S2, it may be possible to reproduce the experi- 
mental spectra. By increasing tq, both transmission co- 
efficient and level density will be affected. It reduces the 
potential barrier, leading to the increase of transmission 
coefficient. Simultaneously, the increase in rg leads to in- 
crease of the leff and vis-a-vis the available phase space. 
We have chosen ro = 1.29 which reproduced the lower 
energy part of the spectra. A similar value was used in 
Ref. [331 . The level density parameter a was taken to 
be A/8 as in previous works [ij, [33].To reproduce the 
present experimental spectra, we have only changed the 
deformability parameters Si and S2, as in the previous 
works [llllllsSl. 



C. Experimental spectra and CASCADE 
predictions 

The measured energy spectra have been compared with 
the respective CASCADE [2^ calculations which have 
been shown in Fig.l. The critical angular momenta, Jcr, 
used in the calculation, are 20, 21, 22, and 23{h) for the 
bombarding energies of 117, 125, 145, and 160 MeV, re- 
spectively [12| . The other input parameters are given in 
Table m The dash-dot-dashed lines represent the results 
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TABLE II. The values difFerent sets of 'deformability' 
parameters: A - obtained from RLDM, and, B - obtained 
by fitting the experimental data (see text). Eiat, E, Jcr and 
Jav are the beam energy, excitation energy, critical angular 
momentum and average angular momentum, respectively. 
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of the CASCADE calculation with the radius parame- 
ter rp = 1.29 and the default values of the deformabil- 
ity parameters, = 3.7x lO""^ and 6^ = 1.1 x 10"^ 
at all beam energies which have been predicted by ro- 
tating liquid drop model [11]. The low-energy parts of 
the experimental a particle spectra are found to match 
with the theoretical spectra, but higher energy parts do 
not. To reproduce the whole spectra, we have followed 
the procedure proposed by Huizenga et al. [sO]- The 
deformability parameters and S2 have been suitably 
optimized to reproduce the experimental spectra, which 
in effect modified the phase space for statistical decay 
by relocation of the yrast line. The results with these 
modified deformability parameters, are shown by solid 
lines in Fig. 1 and the optimized values of deformability 
parameters, 5f and Sf , are given in Table |lTl 

V. DISCUSSION 

It is apparent from Table |lT] that the deformability 
parameters had to be modified from the corresponding 
RLDM values considerably to fit the experimental spec- 
tra. Following the empirical procedure given in Refs. 

[Tsj . it can be shown that the above change in the 
deformability parameters is indicative of substantial en- 
hancement of the deformation of the excited composite 
over the corresponding ground state (RLDM) deforma- 
tion. Before we look further into the origin of such de- 
formation, it may be worthwhile to see if there is any 
other alternative explanation of the observed deviation 
of the experimental spectra from the respective standard 
statistical model prediction. 

Though the present formalism of angular momentum 
dependent level density is largely successful in explaining 
experimental LCP spectra, the magnitude of enhance- 
ment required is quite large and lacking proper expla- 
nation [39I l40j . So an alternative approach, based on 
frozen degrees of freedom, has been proposed, which has 
been shown to reproduce the data quite well [ill, \^ . In 
this approach, it is assumed that the deformation of the 
CN is frozen during the decay, i.e., there is no change 
of shape of the nascent final nucleus; so the phase space 
is calculated using RLDM deformation of the parent nu- 



cleus, rather than that of the usual daughter nucleus. 
This indicates that the dynamical effects such as shape 
relaxation should be taken into account to properly un- 
derstand the phenomenon of particle evaporation from a 
deformed CN. 

The formalism of "frozen deformation" has been ap- 
plied in the present CASCADE calculation. In the case 
of a particle emission, the effect of frozen deformation 
on the energy spectrum may be taken into account in 
the following way. As the a particles are emitted pre- 
dominantly at the initial stage of the decay cascade, the 
deformation may be "frozen" at its value corresponding 
to the highest angular momentum that the CN may have, 
which is w Jcr- So in the present CASCADE calcula- 
tion, the deformation has been kept fixed throughout by 
freezing the value of /e////o, which has been calculated 
using a fixed value of J w Jcr using Eq. [5] with Si and 
62 obtained from RLDM (see Table The resuhs of 
the CASCADE calculation with frozen deformation have 
been shown in Fig. [S] along with the experimental data. 
It is observed that the CASCADE prediction with frozen 
deformation [4lj, |42| for J = Jcr — 23h (solid curve) is 
also in fair agreement with the experimental data. 

To investigate further into the nature of deformation, 
we have computed the values of average effective radius 
parameter, (re//), for all sets of deformability parameters 
(Table |n|. Following [ll|, the effective radius parame- 
ter, re//, is defined as re// = ro\/{l + SiJ'^ + (52 J"'), and 
the average effective radius parameter is calculated using 
the expression, (re//) = '£0'' ''e//(2j + 1)/ Eo'^" (2j -M). 
Then the deviation of (reff) from tq is the measure of de- 
formation, and the deviation of the same from the corre- 
sponding RLDM value may be considered as an enhance- 
ment. The values of re// estimated in the case of 'frozen' 
deformation approach varied between 1.48 (at 117 MeV) 
and 1.58 (at 160 MeV), which are, though higher than 
the corresponding RLDM values ((re//) ^ 1.38 - 

1.42), quite lower than the values obtained with the best 

fitted parameter set {{ref f)^^^*'^'^ ^ 1-52 - 1.66). From 
the above discussion, it may be said that, though the 
present study indicates deformation of the excited com- 
posite which is higher than the corresponding RLDM 
value, uncertainty still remains about the actual magni- 
tude of deformation. This limitation of the present study 
notwithstanding, there is, at least qualitative, indication 
about some enhanced deformation, which may be linked 
with orbiting as one of the contributing factors. 



VI. SUMMARY AND CONCLUSION 

The energy and angular distribution of a particles, 
emitted in the reactions ^^^O (117, 125, 145 and 160 MeV) 
-I- ^^C, have been measured. It has been observed from 
the angular distribution and average velocity plots that 
the a-particles are emitted from fully energy equilibrated 
sources, at all beam energies. The measured energy spec- 
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FIG. 5. Comparison of the experimental energy spectra (cm.) 
of Q-particles (triangles) with the same obtained by theoret- 
ical calculations for beam energy 160 MeV. Solid curve rep- 
resents the CASCADE calculations done using 'frozen' defor- 
mation {lef f /lo) obtained using RLDM deformability param- 
eters for J = 23fi. 



tra have been compared with the same predicted by the 
standard statistical model calculations. It has been found 
that the experimental a particle energy spectra are prop- 
erly explained by CASCADE by optimizing the deforma- 
bility parameters, which are quite different from the re- 
spective RLDM values. This appears to be indicative of 
appreciable change (enhancement) of deformation. How- 



ever, the understanding does not seem to be so simple 
and straightforward when one compares the above with 
the results of CASCADE calculation with frozen defor- 
mation, which has also been shown to be almost equally 
effective in explaining the data. In this case, the effec- 
tive 'frozen' deformation turned out to be smaller than 
that obtained using the optimized parameters. So, the 
uncertainties about the magnitudes of the actual CN de- 
formations notwithstanding, it can, only qualitatively, be 
said that equilibrium orbiting, which is similar to parti- 
cle evaporation in time scale, could also be one of the 
contributing factors for the observed deformation. How- 
ever, the present models are too simplistic to predict the 
actual deformation of the CN; more realistic event-by- 
event Monte Carlo calculations, taking into account the 
initial deformation, spin distribution of the CN and their 
subsequent evolutions, should be performed to have a 
proper understanding of the CN deformation. In addi- 
tion, new experimental inputs (such as measurement of 
deformation from GDR studies [43|) are needed for more 
comprehensive understanding of the process. 
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